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^ ■ Abstract 

. It is shown how quantum field theory at finite temperature can be used to 

QQ . set up self-consistent and gauge invariant equations for cosmological pertur- 

bations sustained by an ultrarelativistic plasma. While in the collisionless 
case, the results are equivalent to those obtained from the Einstein- Vlasov 
J> ■ equations, weak self-interactions in the plasma turn out to require the full 

. machinery of perturbative thermal field theories such as resummation of hard 

. thermal loops. Nevertheless it is still possible to use the same methods that 

I yielded exact solutions in the collisionless case. 

in 
a^ 

I In order to account for the present large-scale structure of the universe such as galax- 

D^' ies, clusters, superclusters, voids, etc., [|[] a cosmological model built on homogeneous and 
isotropic geometries requires certain imperfections in its symmetries. Through the univer- 
sally attractive nature of gravitations, initially small perturbations can grow, in particular 
once that the universe becomes matter dominated and the (maximal) pressure provided 
^ ' by radiation has become inoperative. This picture, which is based on a big-bang scenario, 
• has found dramatic support by the search for and discovery of tiny anisotropies in the cos- 
mic micro-wave background, which in a Friedmann-Robertson- Walker (FRW) universe have 
wavelengths far exceeding the size of the Hubble horizon at the time when this radiation 
decoupled from the primordial matter 0. 

Whatever the origin of these small deviations from homogeneity and isotropy, there 
is a rather long epoch of radiation domination which is thought to be well described by 
a nearly perfect FRW model with metric perturbations evolving in a linear regime. The 
basic equations for these "cosmological perturbations" are nothing else than the perturbed 
Einstein equations. 
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In order to have a close set of equations, these have to be supphed with information 
on the response 5T^" of the energy- momentum tensor to metric perturbations dg""^ . In a 
hydrodynamic approach, this is done by sufficiently restricting the form of ST^" , specifying 
the equations of state, and imposing covariant conversation of the full energy-momentum 
tensor in the perturbed geometry. The simplest case is the one of a perfect (radiation) fluid, 
which has been studied in the pioneering work of Lifshitz P]. Many generalizations have 
since been worked out, and have been cast into a gauge invariant form by Bardeen 0]. A 
modern geometrical justification and generalization has been given recently by Ellis and 
co-authors p. 

A more fundamental description of the behaviour of the primordial matter, which in 
the early universe is mostly a hot plasma of elementary particles, is usually implemented 
through kinetic theory [^]. However, a truly fundamental description eventually has to take 
into account quantum field theory. In the following I shall show that an interesting part of the 
theory of cosmological perturbations can be investigated through the techniques developed 
for quantum field theory at finite temperature [^, namely the case of a weakly interacting 
ultrarelativistic particle plasma. In the limiting case of a collisionless ultrarelativistic plasma 
it turns out to be even possible to obtain exact analytic results where only numerical 
ones where known before; in the case of weak self-interactions one can still find analytic 
results which involve such issues as resummation of hard thermal loops that would be 
very difficult to include in a (quantum) kinetic approach. 

For temperatures T <^ TTipianck it is sufficient to treat the gravitational field as a classi- 
cal background field. The energy-momentum tensor can then be defined by the one-point 
function 

where T[g\ is the effective action functional that contains all the contributions besides the 
classical Einstein-Hilbert action. When derived from this effective action, covariant conser- 
vation of the energy- momentum tensor is automatic and need not be imposed as a constraint. 
The response under perturbations in the metric field is given by 

= / dS'^^Sg-^). (3) 

Hence, ST^^i, is determined by the gravitational polarization tensor (or "thermal graviton 
self-energy" ) 

In particle physics terminology, @ and @ are the sets of one-particle irreducible dia- 
grams with one and two external graviton line(s) in the background field Qf^i, given by the 
cosmological model on which one wants to study the dynamics of cosmological perturbations. 
The concept of thermal equilibrium makes rigorous sense in conformally trivial situations 
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where gfj,u{x) = (T{x)rj^y. This is indeed the case with almost all of the cosmological models 
of interest. If one knows how the effective action transforms under conformal rescalings of 
the metric, then the entire problem of determining the highly nonlocal function (^) (and 
thus the response of the plasma) can be reduced to the evaluation of Feynman diagrams in 
flat space, where momentum-space techniques can be used. In fiat space, temperature can 
be introduced through periodicity in imaginary time, and retarded Green functions in real 
time are obtained by analytic continuation. 

In the high-temperature limit, where all the momenta and masses of the internal particles 
are assumed to be much smaller than temperature, the effective action in fact turns out to 
be invariant under conformal rescalings, so on a curved space with vanishing conformal 
Weyl tensor can be reconstructed by the simple transformation 



= J tl,^,.{k)\^a{y) . (5) 

The Planck mass, which we assumed to be much larger than temperature, does not 
explicitly appear in 11 since we are treating the metric field as classical and no higher loop 
diagrams with graviton self-interactions are involved. So we only need to assume that 
the (zero-temperature) masses of the thermal matter are small compared to temperature 
(i.e., the plasma is ultrarelativistic), and that the relevant momentum scales are likewise so. 
Fortunately, this is just the case of interest with cosmological perturbations. If the latter have 
typical wavelengths of the order of the Hubble horizon, then fc/T ~ V GT^ oc T/mpianck ^ 1- 

One-loop diagrams correspond to coUisionless thermal matter which has only gravita- 
tional interactions. The leading temperature contributions to (^ have been first calcu- 
lated in |jlO| (see also Ref. [|lll]) and turn out to have a universal structure, where only 
the overall factor varies among the various forms of thermal matter according to their en- 
ergy density. It is highly nonlocal and comes with a complicated tensor structure, since 
with T]'^", = 6^, and K'^ = {K^,k), one can build 14 tensors to form a basis for 
W-»l^(K) = p^Ei=i CiiK)Tr'"^{K), see Table |. 

However, H satisfies the Ward identities corresponding to diffeomorphism invariance and 
conformal invariance, and this reduces the number of independent structure functions to 3. 
They can be chosen as 

A{K) = floooo{K)/p, B{K)=flo,%{K)/p, C{K) = fl^,^%K)/ p (6) 

(p = Too), and the Ci..,i4 are determined by the linear combinations given in Table 0. 
The universal result for ultrarelativistic coUisionless thermal matter then reads 

A^^\K) = a;artanhi - ^, B^^^ = -1, C^^^ = 0, (7) 

with UJ = Ko/k. 

Cosmological perturbations can be classified according to their transformation behaviour 
under spatial coordinate transformations [^] as scalar, vector, or tensor, which corresponds 
to compressional, rotational, or radiative perturbations in the plasma. The above 3 inde- 
pendent components of H determine, in certain combinations, the connection between the 
respective perturbations in the energy-momentum tensor and in the metric field. 
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In the radiation-dominated epoch the standard choice is that of a spatially flat Einstein- 
de Sitter model with line element 



ds' = a{T){dT'-di^), a{r) = ^^^T\ (8) 

3 

(po is the energy density when cr = 1), and, given (H), it is moreover natural to decompose all 
perturbations in plane waves, since in linear perturbation theory the different modes evolve 
independently. The problem is thus reduced to a one-dimensional one, and it is convenient 
to introduce a dimensionless time variable 

which measures the (growing) size of the Hubble horizon over the wavelength of a given 
mode (which is constant in comoving coordinates). 

Of all the numerous components of (|I|), only a few are independent by virtue of general 
covariance and turn out to involve only those gauge-invariant combinations of the compo- 
nents of the metric perturbations 5g^y that have been studied by Bardeen |Q . For instance, 
the scalar part of metric perturbations can be parametrized in terms of four scalar functions 



C D, 
Dj A5ij + B^,j 



(10) 



of which always two can be gauged away. Instead of flxing a gauge, we can also use the 
gauge-invariant combinations 

^ = A+''-{D-\b) (11) 
a Z 

li=\B + -B - A) - b --D, (12) 
2 cr a 

where a dot denotes differentiation with respect to the conformal time variable r. 

Each spatial Fourier mode with wave vector k is related to perturbations in the energy 
density and anisotropic pressure according to 

1 1 

5 = ga;^*, T^anis. = ^x^'n. (13) 

Here energy density perturbations 6 are deflned with respect to space-like hypersurfaces 
representing everywhere the local rest frame of the full energy-momentum tensor, whereas 
Tianis. is an unambiguous quantity, since there is no anisotropic pressure in the background. 

Correspondingly, when specifying to scalar perturbations, there are just two independent 
equations contained in (Q). Because of conformal invariance, the trace of (P is particularly 
simple and yields a flnite-order differential equation in x, 

4 12 2 

$" + 1<^' + i$ = f n - -H' (14) 

X 3 3 X ^ ^ 

(a prime denotes differentiation with respect to the dimensionless time variable x). The 
other components, however, involve the nonlocalities of the gravitational polarization tensor. 
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These lead to an integro-differential equation, which upon imposing retarded boundary 
conditions reads 



3)$ + 3a;<l>' = en - 12 f dx' jo{x ~ x'){^' {%') + Y\' {x')) + '^{x - Xq) (15) 



where jQ{x) = sin(x)/x arises as Fourier transform of A{uj) in (0). ip{x — Xq) encodes the 
initial conditions, the simplest choice of which corresponds to ip{x — Xq) oc jo{x — Xq). 

Similar integro-differential equations have been obtained from coupled Einstein- Vlasov 
equations in particular gauges, and the above one can be shown to arise from a gauge- 
invariant reformulation of classical kinetic theory [1^. Usually, these equations were studied 
numerically, but in fact they can be solved analytically . If initial conditions are formulated 
for xo — > 0, a power series ansatz for $ and 11 leads to recursion relations that can be solved 
and lead to an alternating series that converges faster than trigonometric functions. 

This also holds true for the vector and tensor perturbations and when the more realistic 
case of a two-component system of perfect radiation fluid and ultrarelativistic plasma is 
considered [1^ . 




FIG. 1. The energy-density contrast (arbitrary normalization) as a function of x/vr for a colli- 
sionless ultrarelativistic plasma (full line), a scalar plasma with quartic self-interactions X4>'^ and 
A = 1 (dashed line), and a perfect radiation fluid (dotted line). 
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In Fig. 1^, the solution for the energy-density contrast is given in a doubly-logarithmic 
plot (full line) and compared with the perfect-fiuid case (dotted line). In the latter, one 
has growth of the energy- density contrast as long as the wavelength of the perturbation 
exceeds the size of the Hubble radius {x -C 1). After the Hubble horizon has grown such 
as to encompass about one half wavelength [x = vr), further growth of the perturbation 
is stopped by the strong radiation pressure, turning it into an (undamped) acoustic wave 
propagating with the speed of sound in radiation, v = l/VS. The collisionless case is 
similar as concerns the superhorizon-sized perturbations, but after horizon crossing, there 
is strong damping ~ and the phase velocity is about 1. This indeed reproduces the 
findings of the numerical studies of Ref. They can be understood as follows: a energy- 
density perturbation consisting of collisionless particles propagates with the speed of their 
constituents, which in the ultrarelativistic case is the speed of light, and there is collisionless 
damping in the form of directional dispersion. 



While with purely collisionless ultrarelativistic matter, all results are equivalent to 
solving the classical Einstein- Vlasov equations, a quantum-field-theoretical treatment comes 
into its own when self-interactions within the thermal matter are taken into account. In 
a kinetic treatment one could add in a collision term to the coupled Einstein-Boltzmann 
equations, but eventually one would have to abandon the classical concept of a distribution 
function for the thermal matter. A virtue of the above thermal-field-theoretical approach 
is that everything is formulated in purely geometrical terms, without explicit recourse to 
perturbations in the (gauge variant) distribution function. 

In Ref. |]15|, the gravitational polarization tensor has been calculated in a theory 



through order A^^^. The next-to- leading order contributions to HfMuap at order are con- 
tained in the high-temperature limit of two-loop diagrams and their evaluation is straight- 
forward. However, starting at three-loop order, there are infrared divergences which signal a 
breakdown of the convential perturbative series. This is caused by the generation of a ther- 
mal mass oc \/XT for the hot scalars. If this is not resummed into a correspondingly massive 
scalar propagator, repeated insertions of scalar self-energy diagrams in a scalar line produces 
arbitrarily high powers of massless scalar propagators all with the same momentum, and 
thus increasingly singular infrared behaviour (Fig. ^). 





FIG. 2. Two examples of infrared divergent graphs beyond two-loop order. 

However, it is not sufficient to resum this thermal mass for the hot scalars. After all, this 
would break conformal invariance. Indeed, there are also vertex subdiagrams oc AT^ that 
have a similar effect as a self-energy insertion, see. Fig. 03. As in the hard-thermal-loop 
resummation program developed for high-temperature quantum chromodynamics |16|, one 
has to resum also nonlocal vertex contributions. Doing so, the result turns out to satisfy 
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both the diffeomorphism and conformal Ward identities. 

In the low-momentum hmit that is of interest in our apphcation to cosmological per- 
turbations, the function A in (|^ that governs the evolution of scalar perturbations reads 
through order A^/^ 

A = UJ artanh 

UJ 4 



5A 

+ 



+ 



87r2 
87r3 



2 UJ artanh— — uj artanh— 



UJ/ UJ Uj"^ 



3 [uj"^ — 1 — uj\J uj'^ — l) ( UJ artanh— 

UJ 



2 



+6 ( ujVuj'^ — 1 — uj'^ , I UJ artanh— 

UJ Uj"^ UJ /— 2' 

[uj^ — lyi'^ uj^ — I — 1 



(16) 



and similarly complicated expressions arise for B and C, which in the collisionless limit were 
pure numbers. 

The Fourier transform of this expression determines the kernel in the convolution integral 
of ([T5|). At order A^, it can still be expressed in terms of well-known special functions 0, 
whereas at order A^^^ this would involve rather intractable integrals over Lommel functions. 
However, all that is needed for finding analytical solutions is their power series represen- 
tations which are comparatively simple. Given them, it is as easy as before to solve the 
perturbation equations, however one finds that the asymptotic behaviour x ^ 1 eventually 
becomes sensitive to higher and higher loop orders. The reason for this is that higher loop 
orders come with increasingly singular contributions at = ±1 to A{uj), and the large-x 
behaviour is dominated by the latter. This could be cured by a further resummation similar 
to the one introduced for hot quantum chromodynamics in Ref. |T^, but it turns out that 



a particular Pade-approximant based on the perturbative result reflects the effects of this 



further resummation quite well [T^. The result for the density perturbations in a scalar 
plasma with A</)^-interactions and A = 1 are shown in Fig. |l| by the dashed line, where it 
is compared with the collisionless case (full line) and the one of a perfect radiation fluid 
(dotted line). The effects of the self-interactions within the ultrarelativistic plasma become 
important only for x > vr, where the strong collisionless damping is somewhat reduced and 
the phase velocity is smaller than 1. 

A full analysis of scalar, vector, and tensor perturbations in the general case of a two- 
component system containing also a perfect radiation fluid is given in Ref. [|r^. Let us just 
mention one of the more spectacular results, which arise in the case of vector (rotational) 
perturbations. This case has not been investigated much previously, presumably because in 
the perfect-fiuid case there are no regular solutions — rotational perturbations necessarily 
lead to strongly anisotropic initial singularities. This can be explained by the Helmholtz- 



Kelvin circulation theorem ||T9[ which states that in a perfect fluid the circulation around 
a closed curve following the motion of matter is conserved. However, this theorem does 
not apply generally. Indeed, in a (nearly) collisionless medium one can have small initial 
anisotropies in the distribution function that gives rise to a growing vorticity on superhorizon 
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scales W^, which decays after horizon crossing through directional dispersion. In a two- 



component system one can even have such perturbations which do not decay by arranging 
for vorticity in a perfect-fiuid component that is compensated by an initially matching 
one with reversed sign in the nearly collisionless plasma component. Then net vorticity is 
generated by the decay of the vector perturbation in the plasma component, see Fig. |^. This 
is particularly interesting in that vector perturbations generally lead to the generation of 
primordial magnetic fields at the time when the universe changes from radiation to matter 
domination ||20|. Because tiny primordial magnetic fields can act as seed fields for galactic 



dynamos, such rotational perturbations may therefore be of interest with respect to the still 
unsolved problem of the origin of galactic and intergalactic magnetic fields. 
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FIG. 3. Rotational perturbations in a two-component system consisting of 50 % perfect radi- 
ation fluid and 50 % ultrarelativistic scalar plasma (A = 1). Given is the velocity amplitude for 
the total system (full line) and the weakly interacting plasma (dashed line) in arbitary units as a 
function of x/tt. 
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TABLES 



+ (k" + i?^ r?'^'^) + (u^ r/^/^ + u^' 
j.a/3Mi^ = ^i" u/3 + n" + n" n'^ + ul^ u'^ 

j.a/^Mi' = K"" 7?°^ + i?'^ K^' rf" + i?" i?^ + i?'' 77^^ 

r^apnu ^ ^j^p ^oc ^ ^oc ^tx ^ j^u. ^i/^ 

j,apu.v ^ ^ j^fi ^ ^1/ ^ ^, 

= 

rj,a^nu ^ j^f, ^ap ^ j^a ^i^u 

^ ^j^,, ^ j>^, ^_^„a ^ ,,0 ^ J^<^ ^^,„. 



TABLE I. A basis of 14 independent tensors T°''^^'' built from r?'''", = (5^, and 
Ki' = Ki'/k= {uj,k/k). 



ci = iK^A + ^K^B + \C + ^K^ + ^K^ + i 

C2 = + K^B + i^^c + ^^6 + 19^4 + ^^2 _ 1 

C3 = ^2 I f KM + IK^B + i^^c + + 25^4 + ^^2 _ 1 1 

C5 = ^2 I 5^4^ _ 1^2^ _ ic + ^^4 _ 17^2 + ^1 

C7 = U; I - fK^B - IK^C - _ 25^4 _ ^^2 + 1 1 

C9 = (f K2 + f ) + (1^2 + 3) + (1^2 + 1) C- + f ^6 + 9|^4 + |^2 + 1 

Cio = (f + f ) + (1^2 + I) K^S + IK-'C + + + 1^2 + 1 

cn = { (^^2 - i) ^^A - (1^2 + l) S - - - - f } 

C12= te^ + 5^2 + l)^+(5^2 + 2)B + iC+i^4 + iK2 + I 

«i3 = (i^^ + h) K'A - IK'B - iC + l,K^ - - ^ 
C14 = u { fK^A + IK^B +\C- l,K^ + - ^} 



TABLE II. The structure of the conformally covariant gravitational polarization tensor 
m-'^P/p = Y}t^ CiTr""^ in terms of A = floooo/p, B = Uo^,%/p, and C = V'^Vp. 
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